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ABSTRACT 

We present high resolution (R=50,000) spectra at 2.2 /im of 16 young stars 
in the p Ophiuchi dark cloud. Photospheric features are detected in the spectra 
of 11 of these sources, all Class II young stellar objects. The 5 featureless spectra 
consist of 2 Class I, 2 Class 1.5, & 1 Class II. One star, CSS 29, is identified as 
a spectroscopic binary based on radial velocity variations. The radial velocities 
for the remaining sample are consistent with ^^CO and H2CO gas velocities and 
further confirm the membership of the sources in the p Ophiuchi cluster. For 
the 10 spectroscopically single Class II sources, we measure effective tempera- 
tures, continuum veiling, and v sin i rotation from the shapes and strengths of 
atomic photospheric lines by comparing to spectral synthesis models at 2.2 pm. 
We measure surface gravities in 2 stars from the integrated line flux ratio of the 
^^CO line region at 2.3 pm. and the Na I line region at 2.2 pm. Although the 
majority (8/10) of the Class II stars have similar effective temperatures (3530 
K ± 100 K), they exhibit a large spread in bolometric luminosities (factor ~8), 
as derived from near-IR photometry. In the two stars where we have surface 
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gravity measurements from spectroscopy, the photometrically derived luminosi- 
ties are systematically higher than the spectroscopic luminosities. The spread 
in the photometrically derived luminosities in our other sources suggests either 
a large spread in stellar ages, or non-photospheric emission in the J-band since 
anomalous and significant veiling at J has been observed in other T Tauri stars. 
Our spectroscopic luminosities result in older ages on the Hertzsprung-Russell 
diagram than is suggested by photometry at J or K. Most of our sources show 
a larger amount of continuum excess (Fxex) than stellar flux at 2.2 /im (Fk*), 
substantially higher in many cases (rK = Fkcx/Fk* = 0.3-4.5). The derived 
veiling values at K (rx) appear correlated with mid-IR disk luminosity, and with 
Brackctt 7 equivalent width, corrected for veiling. The derived v sin i rotation is 
substantial (12-39 km s^^), but systematically less than the rotation measured 
in Class 1.5 (flat) and Class 1 sources from other studies in Ophiuchus. In four 
stars (Class 1 and 1.5 sources), the absence of any photosphcric lines is likely 
due to large continuum excess and/or rapid rotation if the stars have late-type 
photospheres. 

Subject headings: stars: circumstellar matter — stars: formation — stars: fun- 
damental parameters — stars: Hertzsprung-Russell diagram — stars: low-mass — 
stars: pre-main sequence — stars: rotation-infrared: stars — techniques: spectro- 
scopic 

1. Introduction 

High resolution spectroscopy can be a valuable tool for extracting stellar properties of 
highly extinctcd young stellar objects (YSOs). Spectroscopy offers a more direct method 
than photometry for measuring properties of obscured young stars. Absorption line equiva- 
lent widths and ratios are not dependent on extinction. The intrinsic shapes of photosphcric 
absorption lines, which can be fully resolved in cool stars (T < 5500 K) by high spectral 
resolution (R = A/AA > 10,000-20,000) observations, contain stellar kinematic and tem- 
perature information. Observations of photosphcric lines in obscured stars at high spectral 
resolution permit us to measure (1) precise radial velocities from line shifts, (2) vsini rota- 
tional velocities from line widths, (3) effective temperatures and surface gravities from line 
shapes and line ratios, and (4) the continuum "veiling" by hot dust from the line depths. 

Unfortunately, it is difficult to obtain YSO spectra at optical wavelengths because of 
the heavy extinction in the cores of star forming regions, such as in the p Ophiuchi dark 
cloud (Wilking & Lada 1983). In addition, large continuum excesses from shocks due to 
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accretion of disk material onto the YSO surface often overwhelm the stellar flux in the 
optical and UV (GuUbring ct al. 2000). In the mid-IR, the thermal emission from the 
surrounding disk can also substantially exceed the flux from the photosphere (Fazio et al. 
1976; Lada & Wilking 1984; Bontemps et al. 2001). The standard model of a young star 
surrounded by a circumstellar disk suggests that the photospheric emission should dominate 
the observed emission at near-IR wavelengths. (Strom et al. 1988; Kenyon & Hartmann 
1990). Consequently, the near-IR is the best wavelength region to observe photospheric flux 
from YSOs and can serve as a window through which to observe their physical properties 
and thereby to study earlier stages of star formation and disk accretion than are possible at 
optical wavelengths. 

Clusters are important to study since they represent the dominant mode of star for- 
mation (Miller & Scalo 1978; Shu et al. 1987; Lada & Lada 2003). Since stars form from 
the collapse of dense cores within molecular clouds (Hartmann et al. 1997; Evans 1999), the 
young stars we observe in a cluster have formed from the same material within the parent 
molecular cloud at roughly the same time. Therefore, young clusters present a homogeneous 
sample of stars at similar (though not necessarily equal, Lada & Lada 2003) ages. Embedded 
clusters are of particular interest because they potentially contain the youngest observable 
stars (Wilking & Lada 1983; Lada et al. 1991). 

One of the nearest examples of a young star forming cluster is the dark cloud in the 
p Ophiuchus region (D = 145 pc, dc Zccuw et al. 1999). The compact (d < 1 pc) core of the 
cloud has been well-studied because it contains a large number of embedded low mass stars. 
The first infrared surveys of the Ophiuchus complex resulted in the discovery of dozens of 
embedded YSOs within the more extincted regions of the cloud (Grasdalen et al. 1973; Vrba 
et al. 1975; Ehas 1978a,b). Wilking & Lada (1983), surveyed a 10' x 10.5' region within 
the Ay = 50 boundary (derived from their C^^O contour map) in the p Ophiuchi cloud 
core. With the development of infrared detector arrays, large surveys have been carried out 
to verify cluster membership, detect the presence of disks, and make estimates of the star 
formation efficiency (Greene & Young 1992; Comcron et al. 1993; Strom et al. 1995). Mid- 
IR photometry has been employed to characterize the circumstellar disks around prc-main 
sequence (PMS) stars in this cloud (Greene et al. 1994; Bontemps et al. 2001). A magnitude 
limited (K < 14.0) near-IR imaging survey of a 2.2 pc^ region of p Oph was carried out by 
Barsony et al. (1997). 

In the present study, we will build on the earlier work by using high resolution spectra of 
a sample of PMS stars in p Oph to measure the physical characteristics important to the PMS 
evolutionary stage. High resolution spectroscopy permits us to make direct measurements of 
Teff, log^r, rotation, radial velocity, and excess (non-photospheric) emission. By comparing 
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the effective temperature and gravity or luminosity to evolutionary model tracks on the 
Hertzsprung-Russell (H-R) diagram, we can obtain estimates for YSO masses and ages in a 
sample that is not limited to binary systems. Linking the physical characteristics of young 
stars (singles and binaries, alike) to their PMS evolutionary stage is therefore a primary 
goal. The vsini rotation of PMS stars will give us insights into the relationship between the 
young stars and their circumstellar disks. Quantifying the amount of infrared continuum 
excess allows us to probe the conditions in the disk. Although high resolution spectroscopic 
observations have revealed some of these properties for a few stars in p Ophiuchi (e.g. Bouvier 
et al. 1986; Casah & Eiroa 1996; Greene & Lada 1997, 2000, 2002), no systematic survey 
has yet been completed. 

In § 2, we describe how we obtained our observations and chose our sample in Ophiuchus. 
We present the results of our spectral fitting technique (Doppmann & Jaffe 2003, hereafter 
DJ03) applied to the PHOENIX spectra of our sample of Ophiuchus YSOs, in § 3. In § 
4 and § 5, we discuss the significance of these results, compare luminosities derived from 
photometry and spectroscopy, and put our results in the context of previous work. We 
summarize our conclusions in § 6. 

2. Observations and Data Reduction 

2.1. Selection of the Sample 

The observational goal of this study was to observe a sample of YSOs within the Av = 50 
boundary in the p Ophiuchi cloud core (Wilking & Lada 1983). Wc imposed the following 
selection criteria to determine our sample: (1) We selected only PMS stars that were labeled 
as Class I, 1.5, or II by Luhman & Rieke (1999) (hereafter LR99) or Greene & Lada (1996) 
based on the 2.2 - 10 pm spectral index using photometry of Greene et al. (1994). (2) We 
required a signal-to-noise ratio (S/N) > 20 within one hour to obtain high resolution spectra 
(R ^ 50,000). This resulted in a sample of 18 YSOs with mK < 8.4, listed in Table 1. 

2.2. Spectroscopic Observations 

Prom the sample of 18 young stars in Ophiuchus in Table 1, we obtained near-IR spectra 
of 16 in May 2000 with the PHOENIX spectrograph on the Kitt Peak 4 meter (Hinkle et al. 
1998). Two sources (GY 168 and GY 110) in the magnitude-limited sample were excluded 
bccaTisc there were no nearby optical stars required for guiding on the KPNO 4 meter. The 
positions and alternate source names for the Ophiuchus targets can be found in Barsony et 
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al. (1997). Spectra were imaged on an Aladdin II InSb detector, 1024 x 512 (used) with 
27 /xm pixels. We chose the spectral region at 2.21 nm to measure stellar features in low 
mass stars since this wavelength region has strong photospheric lines in late type stars whose 
shapes and depths are sensitive to changes in temperature (Kleinmann & Hall 1986, DJ03). 

Centered at 2.207 //m, we obtained 0.0095 //m of spectral coverage (the "Na interval" , 
hereafter). With a 1" (4 pixel) slit we achieved a resolving power of R=50,000. 15" nods 
east-west along the slit allowed sequential spectra to be taken in two beam positions on 
the array for removal of foreground sky emission. Targets were kept in the slit by guiding 
on optical stars within the 5' field of the CCD off-axis guider. Single on-chip integrations 
with PHOENIX were 5-25 minutes, resulting in a S/N ~ 40-70 obtained in a single beam 
subtracted pair. We estimated the S/N in our spectra from variations in the photospherically 
featureless continuum shortward of 2.206 /im. Hot stars were observed at the same airmass 
throughout each night for telluric line cancellation. Seven spectral type standards were also 
observed for comparison in our spectroscopic analysis (§3). 

Supplementary high resolution spectra (R ~ 21,000) of the (2-0) ^^CO bandhead (2.2935 
fjLm) of two sources (IRS 2 and GY 314) were generously provided by Tom Greene, fully 
reduced from previously published results (Greene & Lada 1997) using the CSHELL spec- 
trograph (Greene ct al. 1993) on the IRTF 3 meter. The combination of these two data sets 
allowed us to make measurements of surface gravities using the combined spectral windows 
(see § 3.4). 

2.3. Data Reduction 

All data were reduced using IRAF. Source frames at different sht positions were differ- 
enced, and then divided by flat flelds. The flat field source for PHOENIX was a continuum 
lamp that fully illuminated the slit. Bad pixels were removed by interpolation. All spectra 
were optimally extracted using the apall package in IRAF. Telluric absorption lines (H2O 
& CH4) were employed for wavelength calibration using the HITRAN molecular absorption 
database (Rothman et al. 1992), and were removed by division of spectra from hot, early- 
type stars at the same airmass (within 5° of the target) at each slit position. Final spectra 
were obtained by combining the extracted, wavelength-calibrated pairs from the two beam 
positions along the slit. Slopes in the continuum of each source were removed by normal- 
ization to a linear flt across the spectral interval. Figure 1 displays the reduced spectra 
obtained in our Ophiuchus survey. 
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3. Spectroscopic Analysis 

We have developed a spectral fitting technique that provides measurements of physical 
parameters in PMS stars (DJ03). We use information about the line shapes and relative 
strengths to derive effective temperature, v sin i rotational broadening, and continuum veiling 
from fits to spectral synthesis models and are able to measure surface gravities from the 
equivalent width ratios between the Na fines at 2.2 //m and ^^CO lines at 2.3 jim.. Since 
we compare our spectra to photospheric models with different effective temperatures, we 
measure Teg directly without introducing additional uncertainty due to conversion from 
spectral type. Our best fit stellar parameters are determined by the position of the minimum 
in the error space of residuals to the model fits. 

3.1. Effective Temperature, Rotation, and Veiling 

We outline briefly here the technique we have developed for obtaining stellar parameters 
from high resolution near-lR spectra. DJ03 discuss details of the technique, its accuracies, 
and its susceptibility to systematic errors. We use a wavelength region in the K-band at 
2.200-2.210 ^m. (see Fig. 2) to analyze stellar properties by examining the behavior of strong, 
temperature sensitive, photospheric absorption lines (e.g. Na I, Si I, & Sc I). In general, as 
temperature decreases (Teg = 4400-3000 K), the wings of the Na I lines broaden and the Sc 
1 lines grow deeper, while the Si 1 lines decrease in strength. 

For this given wavelength region, we develop a grid of synthetic spectra to which to 
compare our data and determine stellar properties. The synthetic spectra are generated 
from NextGen stellar atmosphere models (Hauschildt et al. 1999) using the spectral synthesis 
program, MOOG (Sneden 1973). The grid of synthesis models covers a range in cool effective 
temperatures {T^s — 3000-5000 K) and PMS surface gravities (logg- = 3.5-5.0). We assume 
solar metallicity in generating the spectra, consistent with the results of Padgett (1996), who 
measures [Fe/H] = +0.08 ± 0.07 for a sample in p Oph (see also the results of Balachandran 
& Carr 1994). 

We add rotational broadening {vsini — 4-40 km s~^ in steps of 1 km s~^) and 2.2 //m 
continuum veiling (rNa = 0-4.0 in steps of 0.1) to these model spectra and compare each 
modified spectrum to the observed spectrum at a given radial velocity. We define the 2.2 /im 
veiling as the ratio of continuum excess at 2.2 /im (the wavelength of the Na interval) to the 
continuum arising from the stellar photosphere (rNa = F2.2^jm excess /F2.2Atm stellar)- A range of 
radial velocities (± 30 km s~^) is tested by shifting the model in increments of 0.8 km s~^ 
(0.5 pixels). We assume the infrared excess is uniform in intensity across our wavelength 
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interval at 2.2 //m. The best model spectrum is selected by finding the minimum RMS of 
the residuals to the fit from selected wavelength sub-intervals that are centered on the Na I 
lines (see vertical dashed lines Fig. 2). 

We start with an assumed gravity of logg' = 3.5, which corresponds to ~l-2 Myr old 
objects in stellar evolutionary models (Baraffe et al. 1998; Siess et al. 2000; Palla & Stabler 

2000), consistent with age estimates of the central embedded cluster in Ophiuchus (Wilking 
et al. 1989; Greene & Meyer 1995; Luhman & Rickc 1999; Bontemps et al. 2001). Next, 
we compare the observed spectrum to models over the entire range of our synthesis grid 
in effective temperature, vsini rotation, and continuum veiling. We then find the best fit 
values of these parameters by examining the RMS deviation of the models from the observed 
spectrum over the regions where lines are present (see D J03) , and then interpolating between 
these models to find the values with the minimum RMS. Figure 2 illustrates our spectral 
fitting technique as it apphes to one of the PMS sources in our sample (GSS 28). 

3.2. Quality of the Model Fits 

The random errors in the model fits are small. We estimate these errors by finding the 
values of the parameters for the minimum RMS when comparing our synthesis grid to target 
model spectra to which we have added Gaussian noise (DJ03). By looking at the variation in 
the values of the parameters at the point where the RMS is minimized for many independent 
Gaussian noise seeds (at S/N = 30 per pixel, Rpix = 240,000), we find (±1 a) uncertainties 
of Teff = 45 K, rNa = 0.04, and vsini = 1.9 km s~^. 

At a S/N of 30, errors in our fits for temperature are dominated by systematic effects. 
When we compare effective temperatures of MK standards derived by our near-IR spectral 
fitting technique to vahies of Tefr derived by converting optical spectral types to temperature, 
we find agreement to a level of ±140 K, la, consistent with the uncertainties involved in 
conversion of spectral type to effective temperature. Tests of the fitting procedure run on MK 
standards to which artificial infrared excesses had been added indicated a small systematic 
overestimate of veiling ArNa= 0.13x(l+ rNa)(DJ03). 

Rotation, temperature, and gravity affect the shapes of photospheric lines in differ- 
ent ways. Fully resolved spectra of transitions with different oscillator strengths in species 
with different energy levels and different abundances can therefore separate the effects of 
changes in vsini, Tes, and logg, even when equivalent width ratios give ambiguous results. 
Nevertheless, there are some combinations of these parameters for which the RMS differ- 
ence between the synthetic spectra in the grid and the target spectrum do not change very 
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quickly. The most notably troublesome combination is temperature and gravity (see § 3.4), 
which we handle by making use of additional information about the (2-0) ^^CO bandhead 
absorption, when available. The quoted random errors were derived holding \ogg fixed at 3.5 
while allowing the remaining three parameters to vary independently and therefore reflect 
the presence of regions in the parameter space with unusually slow variations in the RMS 
difference between observed and model spectra as pairs of parameters vary simultaneously. 

Historically, ultraviolet and optical continuum excesses of T Tauri stars (TTSs) have 
been determined by using the spectra of non-accreting, non-optically veiled stars as a tem- 
plate spectra, as opposed to the synthetic models which wc use as template spectra. In both 
cases the continuum excesses are inferred directly from the diminished line-depths of photo- 
spheric features (e.g. Hartigan et al. 1989, 1991; Basri & Batalha 1990). These veihng studies, 
which use observed template spectra, have demonstrated that the veiling measurements are 
sensitive to the choice of the adopted template. Our method of generating synthetic spectra, 
however, accounts for the lowcr-than-dwarf surface gravity and determines the temperature 
directly from the spectrum, thereby minimizing template mismatches. We demonstrated 
that our synthetic templates are able to reproduce the observed spectra of 9 main-sequence 
stars and 2 luminosity class IV stars (assumed log g = 3.5) without any artificial veiling 
(DJ03). Nevertheless, one may worry that the adopted synthetic model is too simplistic for 
a T Tauri photosphere, which may have star spots and be chromospherically active. How- 
ever, neither of these effects are expected to significantly alter our results since the effects of 
star spots are minimal at near-IR wavelengths (Fernandez & Eiroa 1996; Eiroa et al. 2002), 
where their flux contrast is minimized, and the atomic lines used in the analysis here are 
not chromospherically sensitive, as best we know. Unfortunately we do not have non-veiled, 
non-accreting TTS spectra in hand to directly test synthetic models in these cases. We note 
however, that the optical spectra of non- (or very little) veiled TTSs are reproducible with 
synthetic models (Johns-Krull et al. 1999), which suggests that any possible unaccounted 
for photospheric effects are not likely to bias our inferred properties significantly. 



3.3. Magnetic Fields 

The role of magnetic fields in the evolution of young stars may be significant. Strong 
magnetic fields have been measured in a small number of TTSs (Johns-Krull et al. 1999, 
2001), by modeling the shape of stellar absorption lines that have been broadened by Zccman 
splitting. A typical disk-averaged field strength of 2.5 kG is found from magnetic field 
measurements in six TTSs. 

Lines that are broadened due to Zeeman splitting can introduce a bias into the derived 
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stellar properties if the effects of magnetic fields are ignored in the analysis. Since our spectral 
synthesis models do not have magnetic field strength as a free parameter, we investigate how 
Zeeman broadening changes the derived stellar parameters using synthesis models that do 
account for the effects of magnetic fields on the shapes and depths of the lines in the Na 
interval. We quantify the amount of change in the derived effective temperature, the gravity, 
the vsini broadening, and the veiling that a strong magnetic field can cause by using our 
analysis procedure to derive physical parameters for artificial target stars with a temperature 
of 3600 K, similar to the temperature of most of the sources in our sample. The two targets 
were synthetic spectra provided by Chris Johns-KruU. One synthetic target had no magnetic 
field and the second had an average field of 2.0 kG. We derived Tgs and vsini for these 
targets in the usual way, by finding the position of the minimum RMS difference between 
these artificial stars and our model spectra. The best fit for the target with a 2.0 kG field 
had a lower temperature and a higher vsini value than the target with no magnetic field by 
6% (~1.5 spectral subclasses at MO) and 10% (2 km s~^), respectively. If 2.0 kG is a typical 
disk-averaged magnetic field for the sources in our sample, then our derived temperatures 
are underestimated and rotational velocities are overestimated by a similar factor. 

We also examined how changes in the strengths of lines in the Na interval affect the value 
of loggf derived spectroscopically. We find the best fit to a synthetic target star (see above) 
that has an average magnetic field of 2.0 kG using our models with no assumed magnetic 
fields and allowing only gravity to be a free parameter. We measure a 16% increase in the 
integrated equivalent width across the Na interval when a 2.0 kG field is introduced, which 
causes a decrease in the equivalent width ratio to the magnetically insensitive ^^GO interval 
(see § 3.4). We find that an increase in equivalent width in the Na interval due to the 
presence of a 2.0 kG magnetic field could cause an overestimates in the gravity by ~0.3-0.4 
dex. Since the magnetic field strengths for the stars in our sample are not known, we have 
not corrected for this possible effect. 



3.4. Spectroscopic Determination of Log G 

For late-type stars, the RMS difference in the 2.2 fim region between data and families 
of models in the temperature-gravity plane nmning in the direction from lower temperature 
and lower gravity to higher temperature and higher gravity changes very slowly. If we were 
to adopt logg values that are larger or smaller by 0.5 dex than the value adopted here (see 
§ 3.1), the derived temperatures will change by ± ~200 K and yield similar quality fits see 
(Fig. 6 in DJ03). To resolve this ambiguity, we can use spectra from the (2-0) ^^GO line 
region (2.2925-2.2960 iim). The strength of the ^^GO bandhead is a sensitive luminosity 
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diagnostic in late type stars (Baldwin et al. 1973; Kleinmann & Hall 1986; Ali et al. 1995; 
Ramirez et al. 1997). CO is only mildly sensitive to changes in effective temperature. 

The ratio of photospheric hne strengths between the Na and CO spectral regions depends 
strongly on gravity at a given temperature but is only weakly dependent on temperature 
at a fixed gravity. We can therefore employ this ratio as an indicator of surface gravity in 
the two stars for which wc have measurements of both spectral regions, GY 314 and IRS 
2. We first find the best temperature at an assumed log (7 of 3.5 using our spectral fitting 
routine. For these Tcs values, we then use the ^^CO/Na line strength ratio (Fig. 3) to 
estimate the gravity. We iterate between determinations of Tcfr and logg until the values 
converge (DJ03). This method gives logg values of 4.3 ±0.09 for IRS 2 and 3.5 ±0.10 for 
GY 314. However, the model atmospheres from which the comparison spectra are derived 
are binned in steps of A log 5^=0.5. The best logg^ value for IRS 2 lies in between two models 
on the grid. Since the behavior of the line shapes at 2.2 /im is nonlinear with changes in 
surface gravity, we adopt log 5^=4.0, which has the lowest RMS, rather than interpolating 
between models with log(?=4.0 and log (7=4.5. This choice likely has the effect of lowering 
the determined temperature for IRS 2 by ~180 K, about 1 spectral subclass. 



4. Results 

Figure 1 shows the reduced spectra of all 16 sources we observed in the Na interval. 
The two strongest absorption features are the neutral Na lines at 2.206 and 2.209 /xm. The 
Class I and 1.5 sources (GY 214, GY 6, and GY 252) do not show photospheric features at 
a S/N of ~70. GY 167 (Class II) is also featureless, though at lower S/N. GY 182, reported 
to be early spectral type (B8-A7) by LR99, does not show evidence of photospheric lines at 
2.2 yum, consistent with this early spectral type. 

Ten of 16 sources had detectable photospheric lines (e.g. Na I, Si I, & Sc I) to which 
we could apply our fitting technique. The derived properties are listed in Table 2. These 10 
sources have a radial velocity consistent with that of the p Oph cloud (vlsr of 2.8-4 km s~^, 
Encrenaz 1974; Myers & Ho 1975). The mean radial velocity (vlsr) of these 10 sources is 
3.9 km s~^, with a standard deviation of 1.7 km s~^ about the mean. 

One source, GSS 29, has a radial velocity offset from the cloud of -28.2 km s^^ (see 
Fig. 1). In this case, fits to the photospheric lines are unsatisfactory because the absorp- 
tion lines are too broad and asymmetric, suggesting the presence of a spectrally unresolved 
secondary companion. Subsequent observations of this source have shown it to be a radial 
velocity variable. Although we did not analyze GSS 29 with our spectroscopic fitting routine. 
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we conclude that it is a young spectroscopic binary in p Oph with a mid-K spectral type. 

For eight of 10 sources with fittable hues, we measure cool effective temperatures that are 
within a relatively narrow range (3400-3700 K). The remaining two sources had temperatures 
roughly 1000 K hotter. Spectral types are assigned from the derived temperatures using the 
Sp-Tefif conversion of de Jager & Nieuwenhuijzen (1987). In the 10 sources with good hne 
detections, wc sec a broad range in stellar rotation [vsini — 12-39 km s~^) and veiling at 
K (rNa = 0.3-4.5). 

For two sources (GY 314 and IRS 2), where high resolution spectra including the (2-0) 
^^CO bandhead are available, we employ our full fitting technique (§ 3.4) that uses both ^^CO 
(2.293 yum) and Na I (2.207 //m) line regions to measure the surface gravity, spectroscopically. 
For the remaining eight sources, we assumed a surface gravity of \ogg = 3.5. Table 2 
summarizes all the stellar parameters we measured in our Ophiuchus sample. 



4.1. Luminosities Determined from Photometry 

We derive photospheric luminosities for our sources listed in Table 2 in two ways, both 
based on the near-IR photometry of Barsony et al. (1997). First, we use the published 
J-band apparent magnitudes and assume there is no excess present. We then correct for the 
extinction at J by using the observed J-H colors, adopting the extinction law of Martin & 
Whittet (1990), and applying our measured stellar temperatures to the color-temperature 
relation in Table A5 of Kenyon & Hartmann (1995). Using a distance modulus of 5.81 for 
Ophiuchus (de Zeeuw et al. 1999), we obtain the absolute magnitude at J. We convert to 
bolometric luminosity using the derived effective temperature from our best model fit to the 
data and the bolometric correction factor from Table A5 of Kenyon & Hartmann (1995). 

Our second way of deriving luminosities from the photometry follows the same pro- 
cedure as above except that we use the published K-band apparent magnitudes and re- 
move the fiux that is due to the continuum excess in the K-band that we have measured 
from our best model fits to the spectra. The K-band fiux from the stellar photosphere is 
FK*=FKtot/(l+rNa), where Fxtot is the flux corresponding to the apparent magnitude of the 
source. We then compare the luminosities derived in these two ways by plotting these sources 
against their measured effective temperatures on an H-R diagram together with the model 
tracks of Baraffe et al. (1998) (Fig. 4). 

The dashed lines in Figure 4 connect the stellar luminosities derived from the dereddened 
J-band magnitudes to the luminosities derived from the dereddened K-band magnitudes 
from which we have subtracted the non-photospheric excess using the values of rNa derived 
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from our high resolution spectroscopy. The luminosities based on the J-band photometry 
are systematically higher (by an average of 0.35 in the log) than the luminosities derived 
from the corrected K-band magnitudes. 

4.2. Luminosities Determined from Spectroscopy 

For IRS 2 and GY 314, where we were able to determine loggi directly from spectra (§ 
3.4), we can determine the luminosity independent of the K-band photometry if the stellar 
mass is known. We use the PMS model tracks of Baraffe et al. (1998) to fit our observations 
of log g and T^q to select a unique mass track. Using the relationship between gravity and 
mass, we derive a luminosity: 

L = 47r<jG^ (1) 
9 

Since conversion from gravity to luminosity is dependent on a particular evolutionary 
model, we derive masses and ages for IRS 2 and GY 314 using 4 independent theoretical 
evolutionary tracks (see Table 3). For simplicity, only the Baraffe models are used when 
displaying the luminosities of our sources for which we have spectroscopically measured 
surface gravities. These luminosities along with our measured effective temperatures for 
both sources are plotted on the H-R diagram against the Baraffe stellar evolutionary model 
tracks (Fig. 4). 

5. Discussion 

5.1. Cool Derived Tefr Values and Significant Veiling 

The majority of our sources (8/10) have temperatures near 3500 K, corresponding to 
~M2 spectral type. Seven of 8 of these cooler temperature sources also have temperatures 
derived from lower spectral resolution (R = 800-1200) observations (LR99). In agreement 
with the results presented here (DJW03), LR99 also sec a narrow range in their derived 
temperatures {a = 50 K). The temperatures reported from the lower spectral resolution 
study, however, are systematically higher than the temperatures we derive for these seven 
sources at higher spectral resolution (<Tcfr[LR99] - Tefr[Djwo3]> = 325 K). 

A second systematic difference between our results and LR99 is that the veiling values 
we derive are larger than those determined from low resolution spectra. The median value 
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in the expression: (l+rNa)Djwo3 / (l+rK)LR99 is 1.5 with a range of 1.1 to 2.6. This leads to 
a mean difference in the 2.2 //m continuum veihng of 1.2 from the two independent veihng 
determinations for stars in our sample. Part of the difference between the low and high 
resolution determinations of the 2.2 fim veiling may result from a small systematic error in 
the high resolution rNa fits. When we artificially add continuum veiling to the spectrum of 
an MK standard and then recover this veiling using our fitting routine, the best-fit result is 
systematically larger by ArNa=0.13 x (1+rNa) than the veihng we added to the MK standard 
(DJ03). If a systematic effect of the same size is present in the analysis of Class II sources, 
our analysis of the high resolution spectra overestimates the value of rNa by 0.15-0.7 as rNa 
ranges from 0.3 to 4.5. This effect, however, only accounts for 1/4 of the difference between 
the veiling values determined from our high resolution spectra and the low resolution spectra 
of LR99. Furthermore, tests with higher surface gravity models (log g > 4.5) show less of 
an offset between our recovered veiling values and the continuum excess we added to our 
late-type dwarf standards. 

The largest part of the Tcfr and 2.2 fim veiling difference between our results and those 
of LR99 most likely stems from the differences in spectral resolution. At lower spectral res- 
olution (R < 5,000), where K-band absorption lines from stars with cool photospheres are 
unresolved, measurements of effective temperatures and veiling must rely solely on equiva- 
lent widths. As a consequence, lower resolution measurements cannot distinguish between 
temperature and veiling changes in temperature-sensitive lines because both types of change 
can have similar effects on the equivalent width. The degeneracy between temperature and 
veiling is split when lines are completely resolved and information becomes available from the 
intrinsic shapes of spectral features. We determine more accurate veiling and temperature 
values from the high resolution spectra because we can make use of both the depths of the 
line cores and the shapes of the line wings. 



5.2. Luminosity 

5.2.1. Luminosity and Age Estimates 

The luminosities derived from J magnitudes and from K magnitudes corrected for excess 
emission are plotted against the effective temperatures derived from our high resolution 
spectroscopy in Figure 4. The figure also shows the luminosities determined for IRS 2 and 
GY 314 using our ^^CO/Na hne flux ratio technique (see § 3.4). The spectroscopic technique 
derives logg rather than L. The spectroscopic luminosity therefore depends on the Teff-mass 
relation from the PMS evolutionary models of Baraffe et al. (1998) (also shown in Figure 4). 
The conversion from log g to mass for other PMS models produces comparable luminosities 
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for IRS 2 but results in even lower values for GY 314 (Table 3). The Ophiuchus stars exhibit 
a large spread in luminosity, with 3500 K stars ranging from 0.6 Lq to 5.2 Lq (based on 
the K magnitudes). The Class II sources in our survey correspond to classical TTSs. Since 
these stars are on the convective and primarily vertical portion of their evolutionary tracks 
(Tout, Livio, & Bonnell 1999), this spread in luminosity for stars with similar values of 
Teff may therefore be caused by a spread in age. The photometric results indicate that 
all the stars except IRS 2 he above the 2 Myr isochrone. We note that the young age is 
consistent with values determined for Ophiuchus stars studied spectroscopically and matched 
to different PMS models (e.g. Greene & Meyer 1995, LR99), but caution that interpretation 
of properties at young ages (i.e. < 1 Myr) is subject to considerable uncertainty (Baraffe et 
al. 2002). 

Spectroscopically determined luminosities are largely unaffected by the problems that 
can bias the photometric estimates (e.g. extinction, veiling, unresolved binaries; see below). 
For both IRS 2 and GY 314, the spectroscopically determined luminosities are lower (by a 
factor of 3.2 and 1.4, respectively) than the luminosities determined from the corrected K 
magnitudes. The spectroscopic ages for the two objects also differ sharply (20 Myr vs. 1 
Myr) lending credence to the possibility that a range in ages causes the luminosity spread 
among the 3500 K stars. Although it is difficult to draw conclusions based on only two stars, 
the differences between the photospheric and spectroscopic luminosities and ages lead us to 
make a critical examination of the validity of the luminosities derived from photometry. 



5.2.2. Accuracy of Photometrically Derived Luminosities 

Analyses of the J and K-band photometry both result in higher luminosities than the 
spectroscopic analysis and imply an order of magnitude spread in luminosity for sources 
at 3500 K. While a true spread in source ages may explain the range of luminosities, we 
still need to examine other possible causes for the luminosity spread and also need to look 
for a reason for the offset between the spectroscopic and photometric luminosities and the 
significant offset between the J-band luminosities and the luminosities determined from the 
veiling-corrected K magnitudes. 

1. Near-IR Photometric Variability: Temporal variations in the near-IR energy output 
of a YSO, possibly due to accretion or star spots, will increase the luminosity dispersion for a 
young star sample. However, this effect is usually not large at near-IR wavelengths. Barsony 
et al. (1997) measure J-band photometric variations of Amj > 0.4 (2a) in only 6 of the 79 
sources in their Ophiuchus sample for which they have multiple measurements. Kenyon & 
Hartmann (1995) find a mean variability of ~0.1 magnitudes at K in a sample of PMS stars 



- 15 - 



in Taurus-Auriga. 

2. Unresolved Binaries: Binaries are common in regions of star formation, and probably 
even more common among young stars than field dwarfs (Duquennoy & Mayor 1991; Ghez 
et al. 1993, 1997; Mathieu et al. 2000). Mistaking unresolved binaries for single stars would 
systematically bias the photometry to higher derived luminosities (factor of < 2), introducing 
additional scatter. 

3. Continuum Excess in the J and H-bands: J-band photometry has been commonly 
used to derive stellar luminosity under the assumption that the accretion luminosity is neg- 
ligible at 1.2 fim (Kenyon & Hartmann 1995, LR99). The assumption that one could neglect 
accretion luminosity was tested by Folha & Emerson (1999) who used 1.2 ^m spectroscopy 
to derive J-band continuum veiling in a sample of classical TTSs in Taurus. The average de- 
rived rj value is 0.57, with a maximum of 1.2. The arrow in the upper right corner of Figure 
4 shows that the luminosities derived for our sources in Ophiuchus from J magnitudes could 
agree quite closely with the luminosities derived from the veiling-corrected K magnitudes 
if the average J-band veihng derived by Folha & Emerson (1999) for their Taurus sample 
were removed from the Ophiuchus sources. Unaccounted for non-photospheric emission is a 
plausible explanation for the high luminosities inferred from the J band photometry. 

If veiling is present at H as it is at J and K, and if the veiling flux has an intrinsically 
red color, observers may be overestimating the extinction and therefore overestimating the 
intrinsic flux when they use the observed J-H color to deredden their sources. Meyer et al. 
(1997) examined a population of optically visible classical TTSs. Assuming that the J-band 
excess is 10% of the V band excess derived spectroscopically for these stars, they calculated 
a mean dereddened H-band veiling for their sample of <rH>=0.2. Folha & Emerson (1999) 
use their spectroscopic determinations to derive a mean J-K color of 0.7 magnitudes for the 
veiling flux in their sample. The errors caused by overcorrecting for reddening based on the 
J-H colors would affect our corrected K magnitudes as well as the luminosities derived from 
the J magnitudes. 

5.3. Class I and 1.5 Sources 

We did not detect photospheric lines in any of the Class I or flat spectrum sources 
we observed. The absence of detectable photospheric lines could be due to rapid rotation, 
heavy veiling, diminished line strengths due to higher photospheric temperature or some 
combination of these factors. It is plausible that younger YSOs still undergoing gravitational 
collapse will be rotating rapidly and have significant accretion. Greene & Lada (2000) 



-16- 



conclude that the lack of spectral features in their sample of Ophiuchus Class I YSOs is due 
to heavy continuum veiling (rx > 4.0) and rapid rotation. The broadened lines detected 
in some of their flat spectrum sources is consistent with rapid rotation and is in agreement 
with the vsini = 50 km s~^ measured in one Class I source (YLW 15) which has a veiling 
of rK = 3.0 (Greene & Lada 2002). 

One out of the five sources in our sample that had no photospheric lines detected (GY 
182) was reported to be of early spectral type (B8-A7), while the remaining others had no 
reported spectral type by LR99. It is possible that the other sources have no observable 
K-band lines solely because their photospheres arc hot (Tcfj > 5000 K). In this case, no 
veiling or rotation would be required to produce a featureless K-band spectrum in hot stars. 

To investigate whether rotation and veiling could explain the lack of detectable lines 
in our sample of Class 1 and flat spectrum sources, we set detection limits by finding the 
maximum veiling as a function of v sin i rotation allowed for at least a 3 cr detection of the 
cores of the Na lines (Fig. 5), assuming a temperature of 3600 K. If these stars rotate < 50 
km s~^, then significant veiling (rx > 2.0 for all sources and rx > 4.5 for GY 214) is needed 
in order to hide their photospheric lines below a 3 cr detection limit. 

5.4. ^^CO Line Contributions from the Disk 

We rely on ^^CO equivalent width measurements to determine the surface gravity in 
objects for which we have high resolution observations in two wavelength intervals (2.2 and 
2.3 //m, see § 3.4). While some models predict the presence of ^^CO absorption in the disks 
of TTSs (Calvet et al. 1991), the most definite detections of the CO overtone bandheads 
from disks around TTSs have been in emission (Thompson 1985; Carr et al. 1993; Najita 
et al. 1996, LR99). This ^^CO bandhcad emission is rare and only seems to be associated 
with higher luminosity objects. Wc must avoid using our spectroscopic method on stars 
that exhibit CO emission. The presence of even weak emission may be easy to notice in 
some cases, given that high spectral resolution measurements of WL16 show CO emission 
lines broadened by up to 400 km s^^ (Dent & Geballe 1991). However, Carr & Tokunaga 
(1992), see CO emission from an outflow in the infrared source, SVS 13 with hnes that are 
only 40 km s~^ wide. Vsini broadening in CO absorption hnes of some classical TTSs is 
consistent with stellar rotation rather than Keplerian motion in a circumstellar disk (Casali 
& Eiroa 1996). If the origins of CO emission and absorption at 2.2935 //m are from different 
velocity fields (Keplerian disk vs. stellar photosphere), then at resolving powers of R > 
10,000, the presence of both components would be detectable in the spectrum. Johns-Krull 
& Valenti (2001) find good agreement between their vsini measurements of resolved ^^CO 
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R branch lines in the K-band to the rotational broadening measured in the optical. For the 
two measured sources in our sample, the ^^CO and Na line widths are also consistent. We 
therefore believe that the ^^CO absorption we are using to determine logg accurately reflects 
the line flux arising in the photosphere. 

5.5. Accretion 

Class II YSOs exhibit significant amounts of near-IR emission in excess of the expected 
emission from the stellar photosphere. This excess is apparent in the ncar-IR colors, where 
the colors of dereddened Class II sources form a sequence in the J-H, H-K color plane that is 
distinct from the colors of main sequence and post-main sequence stars (Meyer et al. 1997). 
Measurements of the depths of photosphcric absorption features go farther in showing that 
this near-IR excess is not accompanied by significant line absorption (Greene & Lada 1996; 
Casah & Eiroa 1996; Johns-KruU & Valenti 2001). Greene & Lada (1996) see the near-IR 
excess in Class II sources as part of the overall evolution of the spectral energy distribution. 
They find no detectable fines in the K-band at R = 500 in Class I sources, implying that the 
rising near to mid-IR spectral energy distributions (SEDs) of these objects are accompanied 
by very high near-IR excesses. The Class II objects have significant excesses in the K-band 
while Class III sources have no detectable excess. 

The increasing ratio of excess to photosphcric emission with wavelength makes it clear 
that the near-IR excess, at least at wavelengths beyond 1.2 //m, is not simply a long- 
wavelength extension of the UV excess seen in classical TTSs. The lack of significant line 
absorption from the source of the excess makes it impossible to invoke continuum opacity 
in the gas present in the hotter parts of any accretion disk as the source of the excess 
(Greene & Lada 1996; Calvet, Hartmann, & Strom 1997; Johns-KruU & Valenti 2001). Our 
own data agree with this in showing no absorption lines in the Class I sources and significantly 
weakened but narrow absorption lines in the Class II sources. 

The character of the near-IR excess emission makes it likely that it comes from warm 
dust. The problem with this idea is that it requires both a substantial illuminating lu- 
minosity or viscous dissipation to provide the energy and that some of the dust must be 
extremely hot to produce the shortest wavelength part of this excess. Johns-KruU & Valenti 
(2001) have examined these issues in detail and conclude that, in their sample of classical 
TTSs, the measured accretion rates do not result in sufficient energy from either source for 
any reasonable kind of disk to produce the observed near-IR excess. They argue that the 
strong near-IR excesses require the presence of a significant amount of dust surviving in the 
magnetosphere. 
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For our sample of Class II sources in p Oph, we have measurements of the characteristics 
of the underlying star, as well as a determination of the amount of 2.2 /xm excess from our 
own observations. Measurements of Br7 flux and mid-IR continuum arc also available for 
these sources (LR99, Bontemps et al. 2001). In Figure 6, we compare the 2.2 jim veiling 
(rNa), to various other observed parameters of the Ophiuchus Class II sources. 

The most striking result in Figure 6 is the strong correlation between strength of the 
near-IR veiling (rNa)- and the ratio of disk luminosity calculated from the 7 pm. and 15 
^m. ISOCAM measurements, corrected for the contribution of the stellar photosphere, to 
the stellar luminosity determined from the veiling-corrected K magnitude. The upper left 
panel of Figure 6 displays this relation. The observed correlation shown in the upper left 
panel (correlation coefficient = 0.82) argues that both the near and mid-IR excess depend 
primarily on the accretion rate rather than on variations in the structure of the disk. 

Since most of our sources are faint or invisible at optical wavelengths, it is not possible 
to use the UV excess to directly measure the strength of the emission produced by the 
accretion shock. MuzeroUe, Hartmann, & Calvet (1998) have shown, however, that the Br7 
line emission from classical TTSs can serve as a proxy in determinations of the accretion 
luminosity. The upper right panel of Figure 6 plots a quantity proportional to the ratio of 
Br7 flux to 2.2 ^m. photosphcric flux (the product of the measured Br7 equivalent width and 
l+^Na) against rNa, the ratio of non-photosphcric to photosphcric flux at this wavelength. 
The figure shows a good correlation between EW(Br7) x (l+r^a) and r^a but with a negative 
intercept. One possible interpretation of this result is that both the photosphere and the 
accretion shock heat the dust in the inner disk and magnetosphcrc. At low accretion rates, 
there is still some excess produced by the absorption of the photosphcric emission. As the 
accretion luminosity grows, the amount of reprocessed radiation emerging in the near-IR 
increases while the flux of the photosphere remains constant. 

The two remaining panels in Figure 6 show interesting negative results. The lower 
left panel shows that the veihng is independent of the stellar luminosity. Since 8/10 of 
our sources have a very similar Teff, this result implies that, at least within the Class II 
sources, there is no clear trend relating the veiling to the apparent age. A more reasonable 
interpretation might be that the variation in veiling in a sample with similar ages can swamp 
any age-dependent secular trend in the amount of excess emission. 

The lower right panel tells a similar story. While there is a clear trend of increasing 
near-IR veihng going from Class III to Class I sources, within Class II, there is no signif- 
icant correlation between the veiling and the 2-14 /xm spectral index, the parameter most 
commonly used to distinguish between the SED classes. 
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5.6. Rotation 

The modest rotational velocities of classical TTSs (< vsini > ~15 km s~^, Bouvier 
et al. 1986), despite the strong evidence that these stars were accreting material from cir- 
cumstellar disks, was initially something of a puzzle. Models in which a strong magnetic 
field couples the star to the inner edge of a truncated accretion disk were suggested about a 
decade ago in order to address this problem (Konigl 1991; Shu et al. 1994). In this picture, 
kilogauss scale magnetic fields can lock the stellar rotation rate to the Keplerian rotation of 
the disk at a distance of a few stellar radii where the rotation velocity is modest. Near-IR 
measurements of Zeeman broadening show that some TTSs have magnetic fields of sufficient 
strength for this to work (Johns-KruU et al. 2001). 

One prediction of this magnetic braking picture is that younger stars with high accretion 
rates might cause the disk truncation radius to be smaller than for TTSs and therefore to 
result in a higher rotation rate than for older systems where the magnetic field truncates the 
disk at larger radii. Hartmann (2002) concludes, however, that the timescales for angular 
momentum regulation may be comparable to or longer than the evolutionary timescales. If 
the timescales are indeed comparable, there may be no clear trend in rotation rate with age 
for the youngest objects. 

Greene & Lada (1997, 2002) have argued on the basis of the larger vsini values of a 
few fiat spectrum (Class 1.5) sources, a single Class I source, and the more modest vsini 
measurements of a few Class lis in Ophiuchus, that there in fact is an observable trend 
toward slower rotation as one goes from Class I to Class II. Our additional v sin i data for 
Ophiuchus Class II sources (Figure 7) strengthen this result. The figure plots the near- 
to-mid IR spectral index, the continuous variable used to determine SED class, against the 
V sin i values for our sample, together with the Greene & Lada sources. It shows a clear trend 
toward higher values of v sin i for sources with more positive spectral indices. Among our 
own sample of Class II sources, which have spectral indices between -1.6 and -1.4, however, 
there is little or no correlation with vsini. 

Observers using photometric periods to determine rotation rates have also found evi- 
dence for a relationship between infrared excess (taken in this case as an indicator of accretion 
rate and youth) and rotation. Unlike the IR spectroscopic studies, these optical observations 
examine the rotation behavior for sources with and without near-IR excesses, i.e., those on 
the boundary between Class II and Class III. Edwards et al. (1993) found that most sources 
with periods ~ 8 days have K-band excesses while the majority of sources with 4 day periods 
do not. Choi & Herbst (1996) find that a larger sample of young stars in Orion has a bimodal 
period distribution and argues that the slower peak corresponds to stars with rotation peri- 
ods magnetically locked to their inner disks. The observed bimodal period distribution was 
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later called into question by observations of a larger sample (Stassun et al. 1999), although it 
now appears that the higher mass stars (M > 0.5 M©) do have a bimodal period distribution 
(Herbst et al. 2002). Among the higher mass subgroup, the mean I-K excess of the stars 
with long rotation periods is significantly larger than the excess for stars with short rotation 
periods (Herbst et al. 2002). 

Previous spectroscopic studies indicate that the Class II sources have spun down com- 
pared to the rotation of their presumably younger counterparts (Class 1 and fiat spectrum 
sources). On the other hand, rotation rates based on photometric periods suggest that stars 
without measurable excesses at K (Class 111 sources) have spun up compared to ones with 
significant K excesses (likely to be Class lis). In this context, it is interesting to look at how 
various indicators of disk accretion vary with v sin i for the Class II sources in our sample 
(Figure 8). The upper left panel shows the variation of the ratio of disk luminosity (Bon- 
temps et al. 2001) to photospheric luminosity (this paper) which should reflect the amount 
of heating of the outer disk by photospheric and accretion luminosity. The upper right panel 
shows the variation of the photospheric equivalent width of the Br7 line which also is related 
to the accretion rate (Muzerolle, Hartmann, & Calvet 1998). The lower right panel shows 
the 2.2 nm veiling as a function of vsini. There is no clear trend in any of these panels. 

The absence of a clear relation between vsini and the various accretion indicators in 
Figure 8 can result from either large variations in rotation velocity caused by inclination ef- 
fects combined with random variations in the disk emission characteristics or from a turnover 
from decreasing to increasing rotation rates occurring within the Class 11 population, from 
the comparability of braking timescales and evolutionary timescales, or from the contribution 
of factors other than magnetic locking to the trends in rotation. High resolution near-IR 
spectroscopy cannot produce samples as large as the photometric techniques, but it can give 
us rotation rates for sources across the broadest possible range in SED classes. An expansion 
of our rotation rate sample might still help to clarify the situation regarding the evolution 
of rotation rates. 



6. Conclusions 

We have carried out a high resolution near-IR spectroscopic survey of YSOs in the 
p Ophiuchi dark cloud. We measure effective temperatures, continuum veiling, and vsini 
rotation in ten Class II (T Tauri) sources. For two YSOs where we have high resolution 
spectra at 2.2 ^m and 2.3 fim, we measure surface gravities using a new spectroscopic tech- 
nique we have developed and compare our results to luminosities derived from photometry. 
In five sources where we do not detect photospheric lines, we place limits on rotation rates 
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and the 2.2 //m veiling. 

Prom our measurements, we draw the following conclusions. 

1. With the exception of one source, GSS 29, which has a complex spectrum, all 
the Class 11 sources in our sample have radial velocities within a few km s~^ of the cloud 
velocity determined from millimeter molecular line measurements. This result points out 
the usefulness of high resolution spectroscopy in the near-IR for confirming membership in 
the star-forming cluster. 

2. Most of the brightest Class 11 sources in the core of the p Oph cloud have cool 
derived effective temperatures that fall within a narrow range (< Tegf > = 3530 K ± 100). 
The typical Class II source in our sample has a continuum excess larger than the stellar 
fiux at 2.2 jiTu (< rK > = 2.3). With high resolution spectroscopy we are able to split the 
degeneracy between veiling and temperature which otherwise causes a systematic offset to 
higher temperatures and less veiling, evident when we compare our results to low resolution 
fitting techniques. 

3. Using high resolution spectroscopy of spectral intervals including the Na lines at 2.2 
//m and the (2-0) ^^CO bandhead at 2.3 jim., we derive luminosities for IRS 2 and GY 314 
that are smaller than their photometrically derived luminosities. The ages inferred from 
the spectroscopic measurements are therefore older than the ages derived from photometry. 
Furthermore, these sources do not appear coeval based on the luminosities derived by either 
photometry or spectroscopy. Note, however, that the luminosities may be underestimated 
if the stars have strong magnetic fields. Future derivations of stellar parameters from high 
resolution spetra will need to take the effects of magnetic fields into account. 

4. The vsini rotations of YSOs in p Ophiuchi Class 11 sources (this study), along with 
Class 1.5 and 1 sources (other high resolution studies) indicate that Class 11 YSOs have 
already spun down from earlier evolutionary stages (Class 1.5 & I). 

5. We do not observe any photospheric lines in the Class 1 or 1.5 sources (4 total) we 
observed at S/N > 30, either because they have heavy veiling and/or are rotating rapidly, 
or have hot (Teff > 8000 K) photospheres. 
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Fig. 1. — K-band spectra of Ophiuchus targets observed at high spectral resolution 
(R=50,000) using the PHOENIX spectrometer on the Kitt Peak 4-meter. The spectra 
have been continuum normahzed and are offset for clarity. The first 5 spectra (GY 214, 
GY 6, GY 182, GY 252, and GY 167), show no photospheric features at a S/N of > 30 (see 
§4). 
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Fig. 2. — Spectral synthesis fits to the observed spectrum of GSS 28 (left panels) and residuals 
to each fit (right panels). The variation of line shapes with different effective temperatures 
(top left panel), different vsini rotations (middle left panel), and different veilings (bottom 
left panel) is visible. The best fit to each stellar parameter (holding the other parameters 
fixed) is the middle spectrum in each of the three panels. The best fit model corresponds to 
vsini rotation of 18 km s~^, veiling of 1.2, and effective temperature of 3600 K. The dashed 
lines in the panels indicate the two regions containing photospheric absorption lines where 
the RMS of the residuals is computed in each plot (2.2056-2.2072 //m and 2.2086-2.2092 
urn). 
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Fig. 3. — Isogravity contours as a function of the ratio of ^^CO/Na line flux and effective 
temperature. Equivalent widths were computed from 2.2020-2.2120 fim in the Na interval 
and 2.2925-2.2960 fim in the ^^CO interval. We calculate the surface gravity of 2 Ophiuchus 
sources from their effective temperatures and ratio of photospheric line fluxes measured in two 
regions at high spectral resolution, scaled to the slope of the observed SED at low resolution. 
Teff errors (horizontal error bars) are dominated by systematics in fitting observed spectra 
to the synthesis models (± 140 K, i.e. one spectral subclass). Errors in the measured 
line fiux ratio (vertical error bars) are dominated by the determination of the continuum 
level, a function of S/N. Surface gravities derived using this technique are independent of 
uncertainties in the line-of-sight extinction and continuum veiling that are problematic in 
photometric measurements. 
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Fig. 4. — Differences in derived photometric luminosities of Class II PMS objects in Ophi- 
uchus superposed on theoretical evolutionary tracks (Baraffe et al. 1998). The open squares 
show the luminosities determined from J-band photometry with an assumed excess at J of 
zero. The filled squares show the same sources (connected by dashed hnes) with luminosities 
calculated from K photometry, and accounting for our measured K-band excesses. Addi- 
tionally, the solid circles show luminosities determined from spectroscopy for two sources 
(labeled). The sohd vertical bar in the upper right corresponds to a downward shift in the 
derived luminosity when the J-band veiling increases from zero to rj = 0.57, the average 
veiling value at J seen by Folha & Emerson (1999) in their sample of TTSs. 
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Fig. 5. — Detection limits for the 2.2 /xm Na lines as a function of stellar rotation (vsini) 
and continuum veiling (rx) for the 5 PMS sources that appeared featureless at 2.2 fim. Each 
source would have had detectable Na lines if Tefr=3600 K and the veiling and rotation placed 
it to the lower left of the curve. The minimum rotation and veihng solutions (sohd lines) are 
set by the RMS of the residuals in the Na hne cores to the model fits being 3 a above the 
noise in the continuum. A hard upper limit to v sin i is set by the breakup velocity (dotted 
line). The full range we derive for rotation (vertical dashed lines) and veihng (horizontal 
dashed lines) of the Class II sources (Table 2) with measureable fines in our sample are 
shown for comparison. 
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Fig. 6. — The K-band veilings derived from our spectral survey in Ophiuchus vs. (a) the 
ratio of mid-IR disk luminosity (Bontemps et al. 2001) to stellar luminosity (upper left), (b) 
Brackett 7 equivalent width (LR99) scaled by our derived K-band veiling (upper right), (c) 
stellar luminosity derived from K-band photometry corrected for continuum excess (lower 
left), and (d) the spectral energy index {a = dlog{X F\)/d\) derived from 2 and 14 /xm 
photometry (Barsony ct al. 1997; Bontemps et al. 2001) (lower right). The solid lines in the 
top panels show a linear least squares best fit to the data point. In all panels, the x position 
of the points shows the 2.2 //m veiling (rNa) determined from our spectral fitting technique. 
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Fig. 7. — Ophiuchus YSO comparisons of vsini rotation values vs. spectral index {a = 
dlog{X Fx)/dX) derived from 2 and 14 iim photometry (Barsony et al. 1997; Bontemps et 
al. 2001). Class I and fiat spectrum sources are open symbols. All Class II sources are filled 
symbols. The Class II sources from this study are filled circles. Sources from Greene & Lada 
(1997) and Greene & Lada (2002) are pentagons. The triangles show the measured vsini 
rotation for T Tauri (Class II) sources from optical spectroscopy (Bouvicr ct al. 1986). YSO 
sources with reported vsini measurements from more than one study have points that are 
connected by horizontal dot-dashed hnes. 
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Fig. 8. — Derived v sin i rotation in Class II sources vs. (a) the ratio of mid-IR disk luminos- 
ity (Bontemps at al. 2001) to bolometric luminosity (upper left), (b) Brackett 7 equivalent 
width (LR99) scaled by our derived K band veiling (upper right), (c) stellar luminosity de- 
rived from K-band photometry corrected for continuum excess (lower left), and (d) derived 
continuum excess (veiling) at 2.2 /im from our spectral fitting routine (lower right). 
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Table 1. Ophiuchus Survey Sample 



Source*^ 


SED 
Class 


^ mf 


mn'^ 


mx'^ 


Obs. 
Date (UT) 


approx. S/N 








7 A 




^ /9n /nn 

O 1 zu / uu 




PT 94 


TT 
11 


10 9 


S 9 


u.o 


^ /9n /nn 

o 1 zu / uu 


A^ 


nv 1 fi7 

vj 1 lU 1 


TT 
ii 




o.o 


7 1 
( . 1 


^ /99 /nn 

O / ZZ / uu 


zo 




TT 
ii 


iU.O 


o. I 


7 9 


c, /99 /nn 
o / ZZ / uu 


9'^ 
zo 




TT 
ii 


o.o 


7 7 
I . I 


7 9 
1 .z 


^ /9n /nn 

o / zu / uu 


■^n 

ou 


SR 4 


TT 

11 


Q n 


R n 


7 


1^/91 /nn 

o / zi / uu 


uu 


GY les'^ 


II 


10.3 


8.6 


7.4 






GY 214 


I 


17.2 


12.0 


7.5 


5/20/00 


75 


GY 292 


II 


11.3 


9.3 


7.9 


5/21/00 


60 


GY 110^ 


II 


10.7 


8.9 


8.0 






GSS 28 


II 


9.7 


8.7 


8.1 


5/21/00 


45 


GSS 29^ 


II 


11.1 


9.2 


8.2 


5/21 & 22/00 


55 & 50 


GY6 


I 


13.9 


10.8 


8.3 


5/22/00 


30 


GY 308 


II 


11.5 


9.6 


8.3 


5/21/00 


35 


GY 182 


1.5 


14.1 


10.6 


7.9 


5/20/00 


30 


GY 252 


1.5 


15.2 


11.3 


8.4 


5/20/00 


35 


GY 314 


II 


10.8 


9.3 


8.4 


5/21/00 


40 


IRS 2 


II 


10.5 


9.1 


8.4 


5/22/00 


45 



^Cross reference source names and coordinates are listed Barsony et al. 
(1997) 

^Categorized by LR99 

■^Photometry from Barsony et al. (1997) 

'^Not observed due to lack of nearby guide stars 

^Two observations were made of this source 
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Table 2. Measured stellar parameters of YSOs in p Ophiuchi sample 





-1- eff 


op 


7 I QlTl ? 

U olli I 




vlsr 




T d 


loo" 




(K) 




(km s~^) 




(km s~^) 




(L/Le) 


(cm s~ 


GY 319 


3410 


M3 


26 


0.9 


3.8 


0.14 


1.20 


3.5 


GY 314 


3435 


M3 


22 


1.1 


4.0 


0.83 


0.69 


3.5^ 


GY 23 


3490 


M2 


20 


1.2 


0.1 


1.36 


3.03 


3.5 


GY 292 


3500 


M2 


28 


4.0 


2.8 


1.24 


0.66 


3.5 


EL 24 


3520 


M2 


21 


4.5 


3.8 


1.16 


1.69 


3.5 


SR4 


3585 


Ml. 5 


14 


3.7 


5.7 


0.35 


0.63 


3.5 


GSS 28 


3610 


M1.5 


17 


1.2 


4.2 


0.34 


0.63 


3.5 


GY 308 


3715 


Ml 


17 


1.6 


4.1 


1.15 


0.95 


3.5 


GY 20A 


4460 


K5 


39 


1.6 


6.4 


0.83 


5.19 


3.5 


IRS 2 


4640 


K3 


12 


0.3 


4.2 


0.67 


1.91 


4.3^ 



^Errors are ±140K (DJ03) 

''Based on Tes and adopting the spectral type-Teg relation of de Jager & 
Nieuwenhuijzen (1987) 

'^Determined using J and H photometry of Barsony et al. (1997) our mea- 
sured effective temperatures, the extinction law of Martin & Whittet (1990), 
and the spectral type-color relation from table A5 of Kenyon & Hartmann 
(1995) 

^Determined from the observed K-band flux of Barsony et al. (1997), cor- 
rected for extinction, and the presence of continuum veiling that we derive at 
K 

^Assumed \ogg = 3.5 unless otherwise noted 
^Log g values measured using spectroscopy 
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Table 3. Masses and Ages of p Oph sources with spectroscopically measured 

gravities 



Source 


YSO Evolution 


YSO Mass 


YSO Age 


Luminosity 




Model 


(M/Mq) 


(Myr) 


(L/Lo) 


GY 314 


BCAH^ 


0.5 


1.0 


0.5 


(Teff=3435K) 


SDF^ 


0.34 


1.5 


0.3 


(log g=3.5) 


PS99^ 


0.25 


1.0 


0.24 




DMQS'^ 


0.25 


10 


0.24 


IRS 2 


BCAH^ 


1.0 


20 


0.6 


(Tefr=4640K) 


SDF^ 


1.2 


15 


0.6 


(log g=4.3) 


psgg'^ 


1.0 


20 


0.5 




DMQS'^ 


>0.9 


~70 


>0.5 



^Baraffe, Chabrier, AUard, & Hauschildt (1998) 
''Siess, Dufour, & Forestini (2000) 
^Palla & Stabler (1999) 
^D'Antona & MazziteUi (1997) 



